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Realistic models of amorphous Zr02 are generated in a "melt-and-quench" fashion using ab-initio 
molecular dynamics in a plane-wave pseudopotential formulation of density-functional theory. The 
structural properties of the resulting amorphous models are analyzed, with special attention to 
coordination statistics. The vibrational and dielectric properties of one of these models are then 
investigated from first principles using linear-response methods. The electronic dielectric constant 
and Born effective charges are found to be very similar to those of the crystalline phases. Encour- 
agingly, the predicted total static dielectric constant is ~22, comparable to that of the monoclinic 
phase. This work is motivated by the search for improved gate dielectric materials for sub-0.1 /J,m 
CMOS technology, and may also have implications for Hf02 and for silicates of Zr02 and Hf02. 

PACS numbers: 77.22.-d, 61.43.Bn, 63.50.4-x, 71.23.Cq 



I. INTRODUCTION 

The projections in the International Technology 
Roadmap for Semiconductors call for an effective gate 
dielectric thickness of approximately lnm for 0.1 /xm 
CMOS technology by the year 2006. However, the use 
of current Si02 and silica nitride materials in this regime 
will become problematic because of intolerably high leak- 
age currents. A possible solution, currently under inten- 
sive exploration, is to replace Si02 as the CMOS gate 
dielectric by an oxide having a much higher dielectric 
constant. Use of such a "high-i^T dielectric" would allow 
one to make physically thicker films, hence reducing the 
leakage current, while at the same time maintaining or 
even increasing the gate capacitance. High- AT metal ox- 
ides are at the focus of this effort, with Hf02, Zr02, and 
their mixtures with Si02 showing great promise. 1 ' 2 

In previous work, we studied the structural, electronic, 
and lattice dielectric propterties of crystalline phases 
of Zr02 and Hf02 (Refs. 3-6) using a first-principles 
density-functional approach. We found that the lattice 
dielectric constant depends strongly on crystal phase; for 
example, for Zr02, we find orientationally-averaged total 
e values of 37 and 38 for cubic and tetragonal phases, 
and smaller values of 20, 20, and 19 for the monoclinic 
and for two orthorhombic phases, respectively. We also 
found strong anisotropies, most notably a large in-plane 
susceptibility for the tetragonal phase. Differences in the 
electronic density of states and in band gaps were also 
studied, 5,6 but these are not profound enough to have 
much effect on the purely electronic dielectric constant, 
which remains approximately 5 for all crystal phases. 
These results suggest that novel structural modifications 
of Zr02 or Hf02 might possibly provide a route to en- 
hanced dielectric constants. 

One of the great advantages of using Si02 for CMOS 
technology has been the fact that it forms an amorphous 



oxide, thus allowing it to conform to the substrate with 
enough freedom to eliminate most electrical defects at 
the interface. Materials like Zr02 and Hf02 tend to 
crystallize much more readily than Si02, raising ques- 
tions about whether equally smooth and clean interfaces 
can be formed from such materials. While it is possi- 
ble to form such oxides in an amorphous state by low- 
temperature processes, they tend to recrystallize during 
the thermal treatments needed for later stages of device 
processing. One current avenue of investigation is to con- 
sider crystalline oxides, such as perovskites, that can be 
grown epitaxially on Si. 7 Another is to search for ways to 
raise the recrystallization temperature of Zr02- or Hf02- 
based materials, e.g., by incorporating Si, Al, or N into 
the random network structure. However, very little the- 
oretical guidance is available for such initiatives, since 
almost nothing is known about the structure and atomic- 
scale properties even for pure amorphous Zr02 or Hf02- 

With these motivations, we have embarked on a the- 
oretical study of the structure and properties of amor- 
phous Zr02 (a-Zr02). In this paper, we report the con- 
struction of realistic models of a-Zr02 through "melt- 
and-quench" ab-initio molecular-dynamics (MD) simu- 
lations. This approach, first pioneered by Car and 
Parrinello, 8 combines density-functional theory (DFT) 
with MD into a powerful tool for investigating the physics 
of large systems, especially for liquids and amorphous 
structures where the atomic coordinates cannot be ob- 
tained from diffraction experiments. Using the gener- 
ated models as prototypes for a-Zr02, we investigate the 
atomistic structure, especially the bonding and coordi- 
nation statistics, and then focus in detail on the lat- 
tice dielectric properties. While the calculations have 
been carried out only for a-Zr02, it can be expected that 
many conclusions would hold also for a-Hf02 in view of 
the close chemical and crystallographic similarity of these 
two materials. 

The paper is organized as follows. Section II describes 
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TABLE I. Structural parameters obtained for the three 
Zr02 phases from the present theory, compared with a previ- 
ous theoretical calculation of higher accuracy (see text). Lat- 
tice parameters a, b, and c in A; volume V in A 3 ; monoclinic 
angle /3 in degrees. 



Phase 


Parameter 


Inis Work 


Kel. 3 


monoclinic 


a 


5.098 


5.108 




b 


5.171 


5.169 




c 


5.264 


5.271 




P 


99.49 


99.21 




V 


136.77 


137.40 


tetragonal 


a 


5.037 


5.029 




c 


5.113 


5.101 




d z 


0.041 


0.042 




V 


129.73 


129.03 


cubic 


a 


5.034 


5.037 




V 


127.57 


127.80 



the theoretical approach of our first-principles simula- 
tions, the supercclls used for modeling a-Zr02, and the 
linear-response method used to analyze the lattice di- 
electric properties. In Sec. Ill we present and discuss 
the structural properties of the amorphous models ob- 
tained from the MD simulations, with special attention 
to bonding and coordination. The results of the linear- 
response calculations, such as the Born effective charges, 
phonon modes, and dielectric tensors, are then presented 
and discussed in Sec. IV. Section V concludes the paper. 

II. THEORY 
A. Details of static structural calculations 

In this work, we perform ab-initio constant- 
temperature MD as implemented in the VASP simula- 
tion package, 9 in which the electronic structure is de- 
scribed within DFT 10,11 in the local-density approxima- 
tion (LDA), 12 using a plane- wave basis and ultrasoft 
pseudopotentials. 13 At each molecular-dynamics step, 
the instantaneous Kohn-Sham electronic wavcfunctions 
are obtained using the RMM-DIIS (residual minimiza- 
tion with directive inversion in the iterative subspace) 
method, 14,15 which is found to be particularly efficient for 
diagonalizing the Kohn-Sham Hamiltonian for large sys- 
tems such as the 96-atom supercell considered here. The 
forces at each MD step are calculated as the derivatives 
of the generalized free energy with respect to the ionic 
positions based on the Hcllmann-Fcynman theorem. 16 
The calculated forces are then used to integrate Newton's 
equations of motion for the ionic degrees of freedom via 
the Verlet algorithm. 17 

Because MD simulations on large supercells are quite 
computer-intensive, we have taken some care to make 
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FIG. 1. Relative energies per formula unit versus lattice 
constant a (see Table I) for the monoclinic (□), tetragonal 
(A), and cubic (o) phases of Zr02- 

the calculations as efficient as possible. In particular, we 
have reduced the plane- wave cutoff energy to 15 Ry, and 
the pseudopotential for Zr includes only the outermost 
shells (Ad, 5s) in the valence. To test the adequacy of 
these approximations, we tabulate in Table I the lattice 
parameters calculated with VASP using the current set- 
tings and those determined in Ref. 3, where the cutoff was 
25 Ry and the 4s and 4p Zr levels were included in the 
valence. It can be readily seen that they agree very well. 
In Fig. 1, we show the calculated ground-state energies 
of the three ZrC>2 phases (monoclinic, tetragonal, and 
cubic). The energy ordering (monoclinic lowest, cubic 
highest) is qualitatively and semiquantitativcly correct. 
We thus conclude that the chosen pseudopotentials and 
plane-wave cut-off are adequate to correctly reproduce 
the energetics of the three ZrC>2 phases. 

B. Supercell structures 

In our simulation, a-Zr02 is modeled to have periodic 
boundary conditions with a cubic supercell containing 96 
atoms. This is the same content as a 2x2x2 repetition of 
the monoclinic unit cell, since each monoclinic cell con- 
tains 4 chemical formula units (4 Zr atoms and 8 oxygen 
atoms). In fact, the cubic lattice vectors are only slightly 
distorted from the 2x2x2 monoclinic ones (f3 = 99.21° 
and the lengths of three lattice vectors differ only within 
3%). 3 Thus, we use as a starting point a structure in 
which the lattice vectors are perfectly cubic but the frac- 
tional coordinates are those of the relaxed monoclinic 
structure. (This monoclinic structure remains stable in 
this environment; that is, when it is relaxed under fixed 
cubic lattice vectors, only small further displacements oc- 
cur and the symmetry of the atomic coordinates remains 
monoclinic.) All subsequent MD simulations are done 
using fixed cubic lattice vectors for the supercell, and us- 
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ing single fc-point sampling at the T point of the Brillouin 
zone. 



C. Ab-initio molecular dynamics 

The "melt-and-quench" simulation scheme proceeds 
by carrying out a scries of MD simulations, each at a 
fixed temperature (i.e., in the canonical ensemble). Tem- 
perature T is sequentially increased, step by step, from 
room temperature (T m ) to some hypothetical tempera- 
ture (T max ) high enough to liquefy the system. At each 
step, the structure obtained from the previous MD simu- 
lation is used as the starting point for the next one. The 
process is then reversed to quench the system back to T m . 
For each MD simulation, the system is given enough time 
to reach its thermal equilibrium in order to eliminate its 
correlation to the previous structure. In the case of a- 
ZrC>2, we let the system evolve at temperatures ranging 
from zero to 4000 K, well above the experimental melt- 
ing point of 2980 K. The system is then rapidly quenched 
back to 300 K, followed by a relaxation to zero tempera- 
ture of both atomic coordinates and lattice vectors. 

A canonical ensemble is realized using a Nose thermo- 
stat, which modifies the Newtonian MD by introducing 
an additional degree of freedom such that the total en- 
ergy of the physical system is allowed to fluctuate. 18 ' 19 
In the algorithm of Nose, the extra degree of freedom — 
the Nose mass Q — controls the frequency of the tem- 
perature fluctuations. In principle, the Nose mass should 
be chosen so that the induced temperature fluctuations 
show approximately the same frequencies as the "typical" 
phonon frequencies of the considered system. 20 We know 
from Ref. 3 that the phonon frequencies at the zone cen- 
ter in monoclinic ZrC>2 are in the range of 100- 748 cm -1 , 
corresponding approximately to 40 - 300 fs in period. Our 
experience shows that Q is rather insensitive to the desig- 
nated temperature. In the case of ZrC>2, Q = 0.15-0.20 
has been found to be suitable. For Q in this range, the 
characteristic period of the temperature fluctuations is 
approximately 40 fs, which corresponds the high end of 
the phonon spectrum for monoclinic Zr02- 

It has to be emphasized that VASP currently supports 
only constant-volume MD simulations, i.e., the volume 
of a unit cell is fixed during a "melt-and-quench" simu- 
lation. In order to study the influence of volume on the 
resultant structure, we carry out a series of similar MD 
simulations on supercells of different volumes, as will be 
discussed in details in Sec. Ill A. 



D. Details of the linear-response calculations 

The dielectric properties of the MD-simulated amor- 
phous models of ZrC>2 are calculated by the specialized 
linear-response techniques as implemented in ABINIT, 21 
a simulation package that is also based on DFT with 




1.5 6.0 12.0 18.0 22.5 



Simulation Time (ps) 

FIG. 2. Programmed application of temperature vs. time 
in the "melt-and-quench" simulation. Temperature is in- 
creased from 300 to 4000 K in steps of 1.5 ps, then decreased 
in a reversed sequence, for a total simulation time of 22.5 ps. 

pseudopotentials and plane-wave expansion. ABINIT 
has the capability to calculate response functions such 
as phonon dynamical matrices, dielectric tensors, and 
Born effective charge tensors via a variational approach 
to density-functional perturbation theory 22 ' 23 in which 
phonon displacements and static homogeneous electric 
fields are treated as perturbations. Our linear-response 
calculations are performed within the LDA 12 using the 
Perdew- Wang's parameterization. 24 The Brillouin zone 
is sampled only at the T point. Extended norm- 
conserving pseudopotentials 25 with valence configura- 
tions of Zr(4s,4p,4(i,5s) and 0(2s,2p), and a cutoff en- 
ergy of 35.0 Hartree, are found to provide satisfactory 
convergence. 

III. AMORPHOUS MODELS OF Zr0 2 
A. Generation of amorphous structures 

The "melt-and-quench" MD simulations are schemati- 
cally illustrated in Fig. 2. At each selected temperature, 
the Zr02 system is allowed to run for 1.5 ps to reach 
its thermodynamical equilibrium. The time step for the 
simulations is set to be 3 fs for T — 300- 2500 K, and is 
reduced to 2 fs when T is in the range of 3000 K to 4000 K 
in order to suppress numerical error at higher tempera- 
ture. The total simulation time is 22.5 ps (see Fig. 2). 

As mentioned in Sec. II C, several supercells with dif- 
ferent volume are used in our simulation. The prototypi- 
cal supercell, designated as the p-cell, is adjusted to have 
its volume equal to SVmono, where Vmono is the volume 
of a unit cell of monoclinic Zr02- The side length of the 
p-cell is thus determined as 10.320 a. u.. The other cells, 
designated as the 2-, 4-, 7-, and 11-cell, have their side 
lengths increased by 2%, 4%, 7%, and 11% over that of 
the p-cell respectively. This information is summarized 
in Table II. 

The "melt-and-quench" simulations are performed on 
these five candidate supercells. The results are summa- 
rized in Table III, where Vi is the initial volume of the 
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TABLE II. Supercells used in the simulations. The side 
length a is in atomic units (a.u.). The percentage in the 
second column indicates the increment of a from the p-cell. 



Supercell Type 



Increment (%) 



Side Length (a) 



p-cell 
2-cell 
4-cell 
7-cell 
11-cell 




2 
4 
7 
11 



10.320 
10.527 
10.727 
11.043 
11.467 



supercell that was kept fixed during the MD simulation, 
and Vf denotes the volume resulting from the final zero- 
temperature structural relaxation of coordinates and lat- 
tice vectors at the end of the simulation. The coordi- 
nation number of each atom is determined by counting 
the number of atoms within a cutoff radius of 3.00 A, 
which can be compared with the range of Zr-0 bond 
lengths, 2.035 - 2.233 A, that characterizes monoclinic 
Zr02- 3 It can be readily seen that the volume has a sig- 
nificant effect on the final structure. For the 11-cell, we 
find that the atomistic structure collapses into smaller 
clusters accompanied with "spatial voids" even at room 
temperature (T m ), suggesting that the supercell volume 
is too large to sustain the atomic structure. For the p-, 2- 
and 4-cell, during the "melt-and-quench" simulation the 
coordination numbers suggest some degree of disorder, 
but the final structures recrystallize when the systems 
arc quenched back and relaxed to the ground state. 

Fortunately, a reasonable amorphous model of ZrC>2 
(called "Model I") is realized from the 7-cell. Although 
it has an initial volume of 1346.1 A 3 , a subsequent re- 
laxation to zero temperature reduces the volume to 
1230.8 A 3 , which is approximately 12% bigger than Vi 
of the p-cell. Using Model I as the starting structure, 
we perform a second "melt-and-quench" simulation and 
obtain another model of a-ZrC>2, which is called "Model 
II" hereafter. In the remainder of this section, we will 
mainly focus on analyzing the structural properties of 
two models. 

TABLE III. Summary of structures resulting from the can- 
didate supercells. Vi is the volume at which the MD simula- 
tion is carried out; Vf is the volume after the final structural 
relaxation. Fourth column shows the coordination number 
(CN), or range of CNs, for Zr and O respectively. Last col- 
umn characterizes the resulting structure. 





V (A 3 ) 
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1091.9 
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Crystal 
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Crystal 


7-cell 
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FIG. 3. Distribution of coordination numbers during the 
"melt-and-quench" simulation giving rise to Model I. "Melt" 
and "quench" processes are shown at left and right, respec- 
tively. Simulation temperature is indicated in each panel; Zr 
and O atoms are indicated by solid and dashed bars, respec- 
tively. 



B. Analysis of amorphous structures 

1. Amorphous Model I 

Figure 3 tracks the atomic coordination number (CN) 
of the Zr and O atoms during the formation of Model 
I, starting from the 7-cell in which oxygen atoms have 
CN=3 or 4 and Zr atoms have CN=7. Interestingly, 
the system retreats at T=1500K to a higher-symmetry 
structure with CN=3 and 6 for O and Zr atoms respec- 
tively. A hint of disorder is displayed at T=3000K via 
the incidence of the 5-coordinated Zr atoms and the 2- 
coordinated O atoms. When T is increased to 4000 K, 
the distribution of coordination numbers suggests that 
the system has become strongly disordered, which we 
take as a sign that it has melted. The system is then 
quenched quickly (about 12 ps, see Fig. 2) back to T m . 
From the variation of the coordination numbers, one can 
readily see that the amorphous character has survived to 
T m . A snapshot of this amorphous structure (Model I) 
is shown in Fig. 4, from which we can directly see that 
a reasonable amorphous structure appears to have been 
generated. The numbers on the atoms in Fig. 4 indicate 
the coordination numbers. 

The bond distributions of Zr and O atoms in Model 
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FIG. 4. Structure of amorphous Model I. Zr atoms are 
black; O atoms are white. Coordination number is indicated 
on each atom (for atoms near front or back of view, some 
neighbors may not appear). 
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FIG. 5. Bond length distributions in Model I. Each panel 
shows the number of bonds (frequency) vs. the bond length 
(A) for a particular species ("O" or "Zr") with a particu- 
lar coordination number (e.g., "0 4" indicates 4-coordinated 
oxygen). 



TABLE IV. Bond (Zr-O) statistics in the amorphous mod- 
els I and II. Subscripts of Zr and O denote coordination num- 
bers; N is the number of atoms having the specified coordina- 
tion number. The minimum (L m i n ), maximum (Lmax), and 
average (L) bond length for each type of atom is given in A. 
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Zr 5 


2 


1.90 


2.12 


2.04 










Zr 6 


12 


1.96 


2.32 


2.09 


10 


1.92 


2.66 


2.11 


Zr7 


16 


1.97 
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2.17 
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1.95 
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2.16 


Zr 8 
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2.05 


2.82 


2.25 
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43 
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I are plotted in Fig. 5. While there are 5- and 8- 
coordinated Zr atoms in Model I, it is obvious that the 
6- and 7-coordinated Zr atoms predominate. We can see 
from Fig. 5 that most of the oxygen atoms have a CN of 
3 or 4, except for a few 2- or 5-coordinated oxygen atoms. 

The number of atoms having each CN, together with 
the minimum (L m i n ), maximum (L max ), and average (L) 
bond length for each CN, are reported in Table IV for 
both Model I and Model II. The bond lengths are found 
to be in the range of 2.04 - 2.25 A for L, generally sim- 
ilar to the range of 2.035 - 2.233 A that we found pre- 
viously for monoclinic Zr0 2 . 3 We find that the average 
bond length L decreases monotonically as the coordina- 
tion number increases, consistent with the usual expecta- 
tion, based on chemical principles, of an inverse correla- 
tion between CN and bond strength. The minima (-^min) 
and maxima (L max ) of the bond lengths manifest similar 
behavior. 



2. Amorphous Model II 

As mentioned earlier, a second amorphous structure, 
Model II, is obtained using Model I as the starting input 
structure. Figures 6 and 7 present results for Model II 
in a parallel fashion as for Figs. 3 and 5 for Model I. A 
significant change in going from Model I to Model II is 
the lack of the 5-coordinated Zr and O atoms (Zr5 and 
O5 in Table IV). Interestingly, in Model II, the numbers 
of O3 and O4 atoms have become nearly equal, while in 
Model I there was a substantial difference (iV=43 for 3 , 
18 for O4). The relative populations of Zr 6 and Zr7 do 
not change very much from Model I to Model II, although 
a slight shift from Zr$ to Zr7 has occurred. The average 
lengths of the bonds in the two amorphous models are 
approximately equal. Finally, we point out that the den- 
sity of the system has increased by ^6% in going from 
Model I to Model II. 



5 



□ nn 



FT! r- 

300K 



ML 



2 3 4 5 6 7 8 
Coordination number 



1 2 3 4 5 6 7 
Coordination number 




2.2 2.4 2.6 2.8 
Bond Length 



2.2 2.4 2.6 2.8 3 
Bond Length 



FIG. 7. Bond length distributions in Model II. Each panel 
shows the number of bonds (frequency) vs. the bond length 
(A) for a particular species ("O" or "Zr") with a particular 
coordination number. 



FIG. 6. Distribution of coordination numbers during the 
"melt-and-quench" simulation giving rise to Model II. "Melt" 
and "quench" processes are shown at left and right, respec- 
tively. Simulation temperature is indicated in each panel; Zr 
and O atoms are indicated by solid and dashed bars, respec- 
tively. 

C. Discussion 

We have presented, in this section, two models of amor- 
phous ZrC>2 obtained via a "melt-and-quench" scheme 
using ab-initio MD. Since the volume in our simulations 
is not allowed to vary in any single MD run, certain 
candidate supercells of different volume have been con- 
structed in order to study the effects of volume on the 
final amorphous structures (see Table II). One candidate 
cell (i.e. the 7-cell in Table II) successfully leads to an 
amorphous model of ZrC>2 (Model I). The density of the 
7-cell is about 20% less than the density of monoclinic 
ZrC>2 (pmono), but the final zero-temperature relaxation 
of atomic coordinates and lattice vectors reduces the den- 
sity somewhat, such that the final Model I has a density 
~12% smaller than p m0 no- A second amorphous struc- 
ture, Model II, is generated using Model I as the starting 
structure. The final mass density increases ^6% in going 
from Model I to Model II. 

Although the crystalline phases of ZrC>2 have been ex- 
tensively investigated both experimentally and theoreti- 
cally (Ref. 3 and references therein), there has been very 
little work to characterize the amorphous phase. Exper- 
imental efforts have mainly involved studies of a-ZrC>2 in 
the form of powders 26,27 or thin films. 28 As for a-ZrC>2 



thin films, the resulting structure may be expected to 
depend strongly on deposition and processing conditions 
and on stoichiometric variations and impurities. More- 
over, while ZrC>2 films grown at high temperature are 
typically polycrystalline, there are hints that ZrC>2 films 
grown at low temperature may be a mixture of amor- 
phous and polycrystalline phases. One might also ex- 
pect that the amorphous structure in thin films could 
be quite different from the "bulk" amorphous state that 
our model attempts to simulate. Thus, a direct compar- 
ison of our results with experiments may not presently 
be feasible. In particular, Winterer 28 used a thin-film 
sample (presumably a-Zr0 2 ) with a density of 4.2g/cm 3 , 
approximately 71% of the density of monoclinic ZrC>2 
(Pmono=5.89 g/cm 3 ). The density of his sample corre- 
sponds to that of the 11-cell in our simulation. Accord- 
ing to our calculations, the volume at this density is too 
large to sustain a void-free atomistic structure. On the 
contrary, the densities of two amorphous structures we 
simulated — Model I and II — are about 88% and 94% 
of Pmono- Such a large difference in density would nat- 
urally result in important differences in the local struc- 
ture (e.g., coordination numbers) and the overall physical 
properties. In any case, in view of the potential impor- 
tance of Zr02 and Hf02 in CMOS gate dielectrics, it can 
be hoped that more experimental work in this area may 
soon emerge. 
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IV. DIELECTRIC PROPERTIES OF a-Zr0 2 
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A. Introduction 

The dielectric properties of amorphous Zr0 2 are cal- 
culated using the linear-response features of the ABINIT 
package. Since the structural models generated in 
Sec. Ill B are quite large and have essentially no symme- 
try, the computation of their dielectric properties is quite 
time-consuming. Therefore, we could afford to carry out 
a full study of the dielectric properties of only one of these 
models, and we somewhat arbitrarily chose to focus on 
Model I. 

Before computing the dielectric properties with 
ABINIT, the ground-state structure of Model I obtained 
in VASP is first relaxed again using ABINIT for con- 
sistency. The lattice constant increases slightly from 
10.717A to 10.917A, corresponding to a 5.6% increase 
in volume with respect to the ground-state volume of 

1230.76 A obtained in VASP, while the atomic coordi- 
nates change very little (less than 1%). We use the newly 
relaxed structural coordinates in our subsequent linear- 
response calculations. 

B. Results 

The purely electronic dielectric tensor is computed to 

be 

/ 4.76 -0.03 0.03 \ 
eoo = -0.03 4.62 0.00 . 
V 0.03 0.00 4.54/ 

Clearly this tensor is approximately isotropic and di- 
agonal, as expected for any large supercell containing 
an amorphous material. We obtain an orientationally- 
averaged dielectric constant of 600=4.6. This is only 
slightly smaller than the values obtained previously for 
the various crystalline phases of Zr0 2 (~ 5. 3-5. 7), 3 ' 5 ' 29,30 
confirming our previous conclusion that is fairly in- 
sensitive to the structural phase. 

In order to obtain the lattice contribution to the dielec- 
tric tensor, we first need to compute the phonon mode 
frequencies and effective charges. Beginning with the 
mode frequencies, we present in Fig. 8(a) a histogram 
plot of the phonon DOS. (To be more precise, this is the 
DOS of modes at the Brillouin zone center of the super- 
cell; this accounts for the absence of very low-frequency 
modes in the plot.) As can be seen in Fig. 8(a), the 
phonon modes extend up to about 850 cm" 1 , and ag- 
gregate roughly into three groups with frequencies in the 
range 70 - 300, 300 - 510, and 510 - 850 cm" 1 . The 
overall DOS spectrum does not show the kind of discrete 
features expected for crystalline solids. 

Next, we use the computed mode effective charges Z\ 
and frequencies lu\ of zone-center modes A to construct 
the "infrared activity" shown in Fig. 8(b) and defined 




(a) 



(b) 



150 300 450 600 750 900 

co (cm -1 ) 

FIG. 8. (a) Histogram of phonon density of states (DOS) 
vs. frequency. (Height of each bar is the number of phonon 
modes in the bin divided by the bin width of 10 cm -1 .) (b) 
Histogram of DOS weighted by Z^ 2 /u%. 

as the phonon DOS weighted by Z^/cof. (That is, the 
frequency integral of this function gives the lattice contri- 
bution to the isotropically averaged dielectric constant. 
Here the scalar mode effective charge Z^ a , defined via 
Z^ 2 = J2 a Z\a, is essentially the projection of the atomic 
Z* tensors onto the dynamical matrix eigenvectors as de- 
scribed in Eqs. (6-9) of Ref. 3.) We find that the largest 
contribution to the lattice dielectric response comes from 
modes of low frequency, ^100-250 cm -1 . 

From these ingredients, we calculate the total lattice 
contribution to the dielectric tensor, which is 

/ 20.3 -1.2 0.5 \ 
eiatt = -1.2 17.1 -0.4 . 
V 0.5 -0.4 15.5/ 

Combining with the purely electronic dielectric tensor, 
the total dielectric tensor is then 




25.1 -1.2 0.5\ 
-1.2 21.7 -0.4 
0.5 -0.4 20.0/ 



Once again, we find that these tensors are roughly diago- 
nal and isotropic, which tends to confirm that our struc- 
tural model is indeed amorphous. Taking orientational 
averages, we can summarize our results by observing that 
the purely electronic dielectric constant eoo~4.6 is aug- 
mented by a lattice contribution ei at t— 17.6 for a total 
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FIG. 9. Site-projected phonon density of states g(u>) 
vs. phonon frequency for different atoms types ('02' indicates 
2-fold coordinated oxygen, etc.) in amorphous Model I. 



dielectric constant en— 22. For the crystalline phases, we 
found that eo and ei a tt (in contrast to eoo) are strong func- 
tions of crystal structure, with ei at t ranging from about 
14 to 33, or e ranging from about 19 to 39. 3-6 Here, 
we find that the dielectric constant of a-ZrC>2 (en=22) is 
similar to that calculated previously for the monoclinic 
(eo=20) and two orthorhombic phases (en=20 and 19). 3 ' 5 



C. Decomposition by atom type 

Clearly it is desirable to understand more fully the 
various contributions to the lattice dielectric response of 
the amorphous form of ZrC>2. To this end, we now de- 
compose various lattice properties by "atom type" (that 
is, by chemical species and coordination number) in the 
hope that such an analysis may provide further insight 
into our numerical results. 

We begin by decomposing the total phonon DOS g(u>) 
in Fig. 8(a) into a local DOS on each type of atom corre- 
sponding to the rows of Table IV. The total DOS g(uj) = 
^2\S(u>x — lu) can be decomposed as <?(w) = ^ r g T (uj) 
where the local DOS for atoms of type r is 



]T5>^| 2 <5K-c). 



(1) 



Here uj\ is the eigenfrequency of phonon mode A, r runs 
over atom types (02-05, Zr5-Zr8), and e^ Q is the com- 
ponent of the eigenvector of phonon mode A for atom j 
along Cartesian direction a. The results are plotted in 
Fig. 9. 

We were hopeful that such a decomposition might help 
us understand whether the modes associated with certain 
atom types are systematically much softer or harder. In- 
stead, no clear trends emerge from Fig. 9. (The different 
appearance of the spectra for 05 and Zr8 most probably 
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FIG. 10. (a) Scatterplot of isotropically-averaged atomic 
Z* values (vertical axis) vs. atom type and coordination num- 
ber (horizontal axis) for amorphous Model I. Circles and di- 
amonds denote O and Zr atoms, respectively, (b) Same but 
with 'dielectric activity' (see Eq.(5)) plotted vertically. 



results from the lack of statistics for these atoms, which 
occur in the superccll only once and twice, respectively.) 
A slight reweighting of the spectra towards softer mode 
frequencies for increasing atomic coordination numbers 
can be observed, but this is not a very strong effect. 

Next we test how the atomic Born effective charges 
correlate with atom type. Recall that the Born effec- 
tive charge tensor Z* Q ^ quantifies the macroscopic po- 
larization along Cartesian direction a induced by a dis- 
placement of sublattice i along j3. Due to the amorphous 
nature of our model, the calculated Z* tensors have no 
symmetry, although we find that for most atoms the di- 
agonal elements are dominant. To reduce the large quan- 
tity of Z* data to manageable proportions, we present in 
Fig. 10(a) a scatterplot of just the isotropic averages (i.e., 
one third of the trace) of the atomic Z* tensors sorted 
by atom type for all 96 atoms. 

We find that the Z* values are fairly tightly clustered 
about their average values of 5.0 and —2.5 for Zr and 
O, respectively. These are just slightly smaller than the 
typical values of 5.4 (Zr) and —2.7 (O) computed for the 
crystalline phases, 3,4 but still clearly larger than the nom- 
inal chemical valences of 4 and —2, indicating that signifi- 
cant partial covalent character survives in the amorphous 
phase. 

To investigate the roles of the various atom types fur- 
ther, we have carried out an analysis in which we also 
make an atomic decomposition of the lattice dielectric 
response. 3 To do this, we first decompose the lattice di- 
electric into contributions from pairs of atoms i and j, 
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3? = E g & P) 

where a and (j label the Cartesian directions. Here 




where n\ and e^g are the eigenvalue and eigenvector of 

the force constant matrix <&^ /3 for the phonon mode A, V 

is the volume of unit cell, and R^j — a/3 e j@- We 

then heuristically define the contribution coming from 
atom % to be 

*!S = Y,l ■ (4) 

3 

Finally, we will refer to the trace 

e* = E e - (5) 

a 

as the "dielectric activity" of atom i. 

This quantity is plotted versus coordination number in 
the scatterplot of Fig. 10(b). The results indicate that 
a large number of atoms contribute, but that a notable 
subpopulation of mostly 6-fold Zr atoms are particularly 
strong contributors. 

D. Discussion 

The dielectric constant (eo = 22) of our amorphous 
model is much larger than that of SiC-2 (eo = 3.5) and 
quite comparable to the average dielectric constant of 
the monoclinic phase. Thus, we find that from the point 
of view of this one criterion alone, a-Zr02 is indeed a 
promising high-_ftT dielectric for next-generation gate di- 
electrics. Of course, many other issues need to be ad- 
dressed, not the least of which is the stability (e.g., as 
measured by the recrystallization temperature) of the 
amorphous phase. 

Direct comparison with experiment is difficult, both 
because few experimental measurements on amorphous 
ZrC>2 are available, and because the sample preparation 
procedures may vary and may result in rather different 
amorphous samples. 31,32 We compare here with the re- 
cent work of Koltunski and Devine, 31 focusing mainly 
on the sample deposited at room temperature using no 
rf bias (for which spectra are shown in Figs. 1(c) and 
2(c) of Ref. 31) since the authors judged this sample to 
be their "most amorphous" one. (Upon annealing, the 
samples tended to recrystallize partially to monoclinic 
and tetragonal phases.) These authors reported a dielec- 
tric constant eo (at 100 kHz) of ~ 15 - 18 as measured 
electrically on the as-deposited ZrC>2 thin films. The op- 
tical refractive index measured at 632.8 nm was in the 



range n — 1.83 - 1.85. For this most amorphous sam- 
ple, the absorbance spectrum, obtained using a Nicolet 
Fourier-transform infrared spectrometer, shows some TO 
modes beyond 1000 cm -1 . However, since our calculation 
does not show any phonon modes in this range, we sus- 
pect these may be defect-related features. The sample 
shows broad TO and LO features centered at frequencies 
identified as 410 cm -1 and 693 cm -1 , respectively. 31 (The 
authors also pointed out that use of the Lyddane-Sachs- 
Teller relation would yield an estimate of e ~ 10 for the 
amorphous phase, but this method substantially under- 
estimates eo because the IR measurements did not extend 
to very low frequency.) While our value of 22 certainly 
exceeds the experimentally measured one of 15-18, it is 
roughly the same order of magnitude, and better agree- 
ment might not be expected in view of the fact that the 
nanostructure and stoichiometry of the experimental film 
are imperfectly characterized. 

V. SUMMARY 

We have generated two realistic models of amorphous 
Zr02 via a "melt-and-quench" scheme using ab-initio 
MD simulations. Candidate supercells of different vol- 
ume were constructed in order to study the effects of vol- 
ume on the simulated structures. An amorphous model 
(Model I) was obtained from one candidate structure 
(i.e., the 7-cell), avoiding recrystallization or void forma- 
tion during the "melt-and-quench" sequence. The equi- 
librium density of this model is ~12% less than that of 
the monoclinic phase. The second amorphous structure 
(Model II) was generated from Model I with further pro- 
cessing. We expect both models to be reasonable rep- 
resentatives of a-Zr02- The structural properties and 
the bond-length distributions of Model I and II were 
analyzed, and both models were found to be composed 
mainly of 6- and 7-fold coordinated Zr atoms and 3- and 
4-fold coordinated O atoms. 

The Born effective charges, phonon mode vectors and 
frequencies, and electronic and lattice dielectric tensors 
were then calculated for Model I using linear-response 
methods. The phonon DOS was found to be relatively 
featureless, and the infrared activity spectrum (i.e., the 
DOS weighted by Z^ 2 /uo\) showed a broad peak in the 
range of 100-250 cm -1 in the phonon spectrum. The 
calculations show that the Born effective charges of the 
Zr and O atoms are fairly narrowly distributed about 
5.0 and —2.5, respectively, rather similar to what was 
found previously for the crystalline phases. Both the 
electronic and lattice contributions to the dielectric ten- 
sor were found to be fairly isotropic, as expected for an 
amorphous structure. The calculated electronic contribu- 
tion was about 4.6, slightly less than for the crystalline 
phases. On the other hand, the lattice contribution to 
eo was calculated to be 17.6, rather similar to the aver- 
age value of monoclinic Zr02, though still significantly 
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smaller than for the cubic and tetragonal phases. 

It is difficult to know just how similar our "theoretical 
a-ZrCV is to "experimental a-ZrOV- There are difficul- 
ties on both sides of the comparison. On the theoretical 
side, several approximations have been introduced. First, 
the modest size of our model supercell (2x2x2) is a se- 
rious approximation, even though it is essentially at the 
limit of current computational capabilities. However, the 
nearly isotropic form of the lattice (ei a tt) and electronic 
(too) dielectric tensors tends to confirm that we have lost 
the memory of the initial crystalline structure and ob- 
tained an amorphous one, in spite of the supercell size 
limitation. Second, recall that the 4s and 4p electrons 
were included in the valence in the Zr pseudopotential 
for the MD simulations, but not for the linear-response 
calculations. This additional approximation could affect 
the accuracy of the computed dielectric properties of our 
amorphous structure. However, we believe this approxi- 
mation is not very significant since the structural calcu- 
lations on the monoclinic phase of ZrC>2, which is rather 
"disordered" itself, give nearly identical results using the 
two Zr pseudopotentials. Experimentally, it is possible 
that the amorphous state manifests itself differently in 
bulk and thin-film forms. Indeed, because ZrC>2 is es- 
sentially a poor glass-former, it can only be prepared 
in amorphous form by low-temperature deposition or in 
other non-equilibrium conditions, and it is therefore to 
be expected that the properties of such samples may 
vary strongly with preparation conditions. Motivated 
by applications, much work is now under way to explore 
whether chemical substitution (e.g., by addition of Si, Al, 
or N) may help stabilize the amorphous phase and raise 
its recrystallization temperature, and ultimately it will 
be important to extend the theory to such substitutions 
in the future. Nevertheless, we believe that the present 
computational investigation can serve as an important 
first step in answering some of the many open questions 
about the structural and chemical properties of this class 
of materials. 
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